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Abstract. We present a systematic study of capillary filling for a binary fluid by using a mesoscopic lattice 
Boltzmann model for immiscible fluids describing a diffusive interface moving at a given contact angle with 
respect to the walls. The phenomenological way to impose a given contact angle is analysed. Particular 
attention is given to the case of complete wetting, that is contact angle equal to zero. Numerical results 
yield quantitative agreement with the theoretical Washburn's law, provided that the correct ratio of the 
dynamic viscosities between the two fluids is used. Finally, the presence of precursor films is experienced 
and it is shown that these films advance in time with a square-root law but with a different prefactor with 
respect to the bulk interface. 

PACS. 47.ll.-j , - 68.08.Bc , - 68.15.-he , - 81.15.-z 

1 Introduction put forward since decades [HH]. However, since capillary 

filling is a typical "contact line" problem, where subtle 



Capillary filling is an inrportant subject of research for 
its relevance to microphysics and nanophysics [l][2] , which 
deal with the development of nrodern nanofluidic devises, 
as well as with various applications of porous materials. 
The fundanrental physics of capillary filling has been widely 
studied from the pioneering works of Washburn [3] and Lu- 
cas [1] and the basic equations of such problem have been 



non-hydrodynamic effects take place at the contact point 
between liquid-gas and solid phase, it remains an impor- 
tant subject of experimental and theoretical present re- 
search at micro and nano-scales [71I51[^[TU] . Usually, only 
the late asymptotic stage is studied, leading to the well- 
known Washburn's law, which predicts the following rela- 
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2 S. Chibbaro: Capillary filling with pseudo-potential binary Lattice-Boltzmann model 

tion for the position of the interface inside the capillary: apphcations which can not be described within a pure 
2/,\ 2(n\ _ "fHcos(9) continuum framework but do not require a fully atomistic 



3fi 
where 7 is the surface tension between liquid and gas, 9 is 
the static contact angle, /i is the liquid dynamic viscosity, 
H is the channel height and the factor 3 depends on the 
geometry of the channel (here a two dimensional geometry 
given by two infinite parallel plates separated by a distance In this work, we propose to use the Lattice-Boltzmann 



treatment. Many interesting and challenging phenomena 
can be placed in this category, like impact of a droplet 
on solid walls, capillary dynamics and water repellency on 
structured surfaces ^0112 11122] , 



iJ - see fig. [T|). model (LBM) with a pseudo-potential interaction, first 

Experiments of such phenomena are quite difficult, proposed by Shan and Chen [23ill3j and already success- 

mainly at nanoscopic scales, and thus numerical Simula- fully used in many applications in micro/nanoscience |241 

tions are an essential tool for the investigation of this prob- [^[^ . Specifically, we use the model proposed for the de- 

lem. One possible approach is based upon an atonfistic scription of immiscible multi-fliuds, limiting ourselves to 

description, where Newton's equations are integrated for a binary mixture. For this kind of model, some numeri- 

a set of molecules, which interact typically via a Lennard- cal simulations have been worked out, although in sim- 

Jones potential [11. This idea is behind the molecular plified configurations |28|. Nevertheless, this kind of mod- 

dynamics (MD) approach, which is regarded to reproduce els seems to assure a good compromise between efficiency 

rather well the reality. However, its computational effi- and accuracy. On the one hand, these models are hydro- 

ciency is very low and macroscopic scales are very hard dynamic and, therefore, assure the correct macroscopic 

to reach. Therefore, different models have to be inves- model, since there are many evidences that the capillary 

tigated. The Lattice-Boltzmann (LB) method p!2l[T4lfT5] imbibition remains hydrodynamic down to nanoscales [29l 

has been materialising during the last decades as an ef- [SU]. On the other hand, the formulation in terms of a dis- 

ficient and powerful mesoscopic way to simulate complex cretised Boltzmann equation allows a high efficient and 

flows. Recently, this approach has been also assessed for easily parallel algorithm and the possibility of introduc- 

the description of fluid-solid interactions, which have a ing in a direct manner some ingredients which are taken 

major role for micro- and nanoscopic applications [161117] . from more microscopic approaches. The study of capil- 

In particular, different models for slip boundary conditions lary imbibition was recently attempted via Shan-Chen 

and wetting properties have been put forward [3Tl[T8ll24l one-component (multi-phase) LB method [27], which has 

[TM77] . Given that the LB phcnomcnological description been found to reproduce correctly the basic physics of the 

is not based on atomistic details but on average prop- capillary flUing. Nevertheless, this was found to be true 

erties, this mesoscopic model looks promising for those only provided that the evaporation-condensation effect is 
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negligible, that is when the gas phase density is neghgi- At a time step), and f^ {x,t) the corresponding equi- 

ble with respect to the hquid one 27J. In this sense, the librium function. The coUision-tinie is related to kine- 

binary LBM seems to present a more general framework matic viscosity by the formula v^ = h{Tk — ^)- For a two- 

with respect to the multi-phase one, since the evaporation- dimensional 9-speed LB model (D2Q9) f^ {x,t) takes 

condensation effect is very low for this model. the following form [15]: 



In this paper, we show the first evidence that the LBM „ 
chosen can satisfactorily reproduce the physics of the cap- 
illary filling in realistic configurations (for example Water- 
air mixtures) with a proper choice of the physical parame- +2*^'^* fe' g'^fc fc — ... 

, , 1 fk(ei) _ (1 ^ Oik)nk 1 eq /r\ 

ters. A given contact angle can be easily imposed through Ji — ^tj r —n^Ci ■ Uj, (O) 

the introduction of phenomenological force between fluids -| — Uhici ■ u?'')^ nkuf^ ■ uf' for i=5. . .8 



fo = OLkUk - ^J^fcU^."^ • u^' (3) 

It'-'' - ^^^^ + ^^^.c. • u- (4) 

' 1'^'?^2 eg eq 



and walls. The interesting situation of a complete wetting 
is also supported by this approach. The expected Wash- 
burn's law is found for the displacement of the interface 
between the two fluids. Moreover, precursors films appear 
to be present and to follow a Washburn's law, with a dif- 
ferent pre-factor. All the methodological features are de- 
tailed. 



24 

In the above equations c^'s are discrete velocities, defined 
as follows 

0,i = 0, 



Ci = < 



(j:-l)7r . (i-l)7r\ . -I , 



V2 [cos[^-^ + f ], szn[i^ + |] j , i ^ 5 _ 8 

(6) 

in the above, ak is a free parameter related to the sound 
speed of the kth component, according to (c^)^ — 1(1 ~ 
ctk)', nk — X^i fi is ths number density of the fcth com- 
In this work we use the muhicomponent LB model pro- ponent. The mass density is defined as pk = m^Tifc, and 



2 Model 



posed by Shan and Chen [23 . This model allows for dis- ^^^ A^id velocity of the kth fluid u^ is defined through 
tribution functions of an arbitrary number of components, Pk^k = "^fcZ^iCi/j , where m^ is the molecular mass of 
with different molecular mass: t^^ ^t^i component. The equilibrium velocity u^' is deter- 

mined by the relation 



ftix + c,At,t + At)-ft{x,t)^ (2) 



At r 

Tk 



f^(x,t)-,f^^^''\x,t) 



Pfeuf = PkV + TkFk (7) 



where f^{x,t) is the kinetic probability density function where U is the common velocity of the two components, 
associated with a mesoscopic velocity C; for the fcth fluid. To conserve momentum at each collision in the absence of 
Tk is a mean collision time of the fcth component (with interaction (i.e. in the case of Ffe = 0) U has to satisfy 
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the relation map of the value of the static contact angle versus the 

(v^ Oi ut \ / v^ Di \ value of the interaction strength g^ has to be obtained. 

I / \ I / To this aim, we have carried out several simulations of 

The interaction force between particles is the sum of a , ,• i i . .. i i ; ^^ r i-rr , i r 

a static droplet attached to a wall tor dinerent values ot 

bulk and a wall components. The bulk force is given by i „^ ii „„i -, , . , . , , , , r , i 

gw |3iy38j . in particular, m our work, the value ot the 

s 

Fn^^(x) = —^dx)'S^^S^Gi,i,^i,{x')(x' — x) (9) static contact angle has been computed directly as the 

slope of the contours of near-wall density field, and inde- 



Gfcfc — ^ 



where Gj.j, is symmetric and ^k is a function of rik ■ In our 
model, the interaction-matrix is given by 

9k-kAx' -x\ = 1, 

g,,/4, \x' -x\ = V2, (10) 

0, otherwise. 

where g^.^. is the strength of the interparticle potential 
between components k and k. In this study, the effective 
number density ^k{nk) is taken simply as ^k{nk) = nk- 
Other choices would lead to a different equation of state 
(see below). 

At the fluid/solid interface, the wall is regarded as a 
phase with constant number density. The interaction force 
between the fluid and wall is described as 



pendently through the Laplace's law, AP = 27ras0 rpj^^ 
value so obtained is computed within an error ^ 2% — 3%. 
Recently, a different approach has been proposed, which is 
able to give an "a priori" estimate of the static contact an- 
gle from the phcnomcnological parameter [33] • Neverthe- 
less, we have preferred to retain our "a posteriori" method 
for its simplicity and etflciency. 

In a region of pure /cth component, the pressure is given 
by pk — {c^)^mknk, where (c^)^ = |(1 — a^). To simulate 
a multiple component fluid with different densities, we let 
(Cg)^r?T,fe = Cq, where Cq = 1/3. Then, the pressure of the 
whole fluid is given by p = c^ J2k ^k + § Efe fe dk^k'^k^'k 



which represents a non-ideal gas law. 
F2fc(x) = -nk{x)'^gkwPwix'){x' - x) (11) The Chapman-Enskog expansion [15] shows that the 

x' 

fluid mixture follows the Navier-Stokes equations for a sin- 
gle fluid: 



dtp + V-{pn)^0, (12) 

p[dtu + [u ■ V)m] = -VP + F + V • {p{Vu + uV) 



where p^ is the number density of the wall and g^^ is the 

interaction strength between component k and the wall. 

By adjusting g^^ and pw , different wettabilities can be 

obtained. This approach allows the dcflnition of a static 

contact angle 9, by introducing a suitable value for the wall 

density p^ [31], which can span the range 6* e [0° : 180°]. where p — ^j, pk is the total density of the fluid mixture 

It is worth noting that, with this method, it is not possi- the whole fluid velocity u is defined by p\i — ^^, pkV-k 4 

ble to know "a priori" the value of the contact angle from ^ ^j, F^ and the dynamic viscosity is given hy p, = pv = 

the phcnomcnological parameters. Thus, an "a posteriori" ^j, p^ — '^kiPk^k)- 
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2.1 Theoretical analysis nel, namely: 

u 

'^yy> ^ ^jpyyH - y) 

For the geometry set up considered in this work, see fie. , _ rii ^ i tt r i r \ ■ i ■ , r ^ 

where u = ^ — 1/J^ J ayu[y) is the mteriace mean veloc- 

[U we propose a general model which can be integrated 

ityj 

(in this paper we have used the software "Mathematica" ) o i « "" 

H 

to provide a semi-analytical solution to be compared with 



dyU\y^Q = 6— (14) 



Inserting this expression in the viscous force and following 
the notation of fig. ([1]), one obtains the final expression (see 
also [32l[27l for a similar derivation: 



numerical results. 

The Washburn's law ([1]) holds when inertial forces can 
be neglected with respect to the viscous and capillary ones 

[5]. This is not true in the early stage of the filling pro- {Pg{L - z) + piz)z + {pi - Pg){z) 
cess, typically a few nanoseconds for micro-devices, where 2— -^ [(Mq(^ ^ ^) + ^^lz)] (15) 



2 



strong acceleration drives the interface inside the capillary. 
Another important effect is the "resistance" of the gas oc- 
cupying the capillary during the liquid invasion. This re- 
sistance is due to the finite value of the gas density and is a 
sensitive issue for numerical simulations, because reaching 
the typical 1 : 1000 density ratio between liquid and gas 
of experimental set up, represents a challenge for most 
numerical methods. In order to take into account both 
effects, inertia and gas dynamics, one may write down 
the balance between the total momentum change inside 
the capillary and the force (per unit width) acting on the 
liquid+gas system: 



dizMit)) , , 

^ cap 1^ -'- VIS \-^^) 



In the above equation for the front dynamics, the terms in 
the LHS take into accoiuit for the fiuid inertia. Since these 
terms are proportional either to the acceleration or to the 
squared velocity, they become negligible for long times, 
as the velocity becomes negligible too. The right hand 
side takes into account for the forces. The presence of two 
viscous forces related to both liquid and gas components 
indicates that a gas resistance is felt for Pg/ p-i ^ 0. Thus, 
Washburn's law ^ is correctly recovered asymptotically, 
for t -^ oo and in the limit Pg/ pi — > 0. 

3 Numerical Results 



dt ^ 

In figure [1] the set-up of our numerical experiments is 

where M{t) ~ Mg + Mi is the total mass of liquid and gas sketched. In order to simulate a realistic capillary we im- 

inside the capillary at any given time. The two forces in the pose that the bottom and top surfaces are coated only 

right hand side correspond to the capillary force, Fcap = in the right half of the channel with a boundary con- 

2"fcos{9), and to the viscous drag Fyis = —2(pg(L ~ z) + dition imposing a given static contact angle [31J ; in the 

piz)dnu{0). We consider the Poiseuille solution for a chan- left half we impose periodic boundary conditions at top 
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and bottom surfaces in order to mimic an "infinite reser- constant in time, a simple analytical solution of equation 

voir" . Periodic boundary conditions are also imposed at (flS)) can be obtained: 

the two lateral sides such as to ensure total mass conser- 

vation mside the system. At the solid surface, bounce back; QfiL 

boundary conditions for the particle distributions were im- 

where zq is the starting point of the interface at the be- 

posed. It is important to note that, in the following, we 2 

ginning of the simulation, td = j^ is a typical transient 

shall indicate the first component with an index 1, and 

time. The analytical solution shows an asymptotic linear 

the second with 2, for the sake of clarity, while in figure 

dependence in time for the front displacement z{t). This 

[1] they are respectively indicated by 1 (liquid) and by g 

behaviour is correctly reproduced by our LB model, as 

(gas). The dimension of the channel have been chosen to 

shown in figure [3l where the numerical solution is com- 
be: L = 200 Az ; H = 40 Az. It is important to note 

pared against the analytical solution given by eq. (|16p . 

that in LB methods the interface S^ is diffuse {d(, sa 4 — 5 

Recently, the issue of simulating components of dif- 
grid points in present study) and it has been pointed out 

ferent densities has been addressed for problems of binary 
that meaningful results can be obtained only when the ra- 

diffusion [3S]. However, in our knowledge, there is no com- 
tio S^/H is large enough, S^/H -- 10 [27]. Therefore, the 

plete study of a dynamical phenomenon with components 
size of the channel has been chosen such as to fulfill this 

with different densities. Following a recent work [36], it is 
condition. 

pointed out that multicomponent LBE approach describes 

The first issue to address in order to have a realistic two isothermal perfect gases, in mutual interaction, which 

simulation of a capillary filling, is the correct description of conserve their isothermal thermodynamics. In an isother- 

the physical properties of the system. It is well known that nial gas, the sound speed is inversely proportional to the 

in LB methods it is difficult to maintain realistic density molecular mass ci = -^^^, where T is the gas adiabatic 

ratio between different species. This is particularly true for constant and R the universal gas constant. In our model, 

binary LB methods, where almost all previous simulations molecular mass can be changed for each component to- 

have been limited to density ratio equal to 1 [^51[5T||331[3^ ■ gether with the value of the parameter ak in order to 

As a first vahdation, we present the results obtained for change the value of the sound speed in eqs. ©-([S]). For 

a dynamic configuration where fluid and gas have equal a binary mixture, the two sound speeds are related as 

density (pi = p2 = 1) and equal viscosities (^1 = /i2 == c] ^ JW'^l- This relation allows to fix the value of 02, 

^ = i). The contact angle is chosen to be 6* « 40, that once ai has been chosen through the relation 
is giw = 0; g2w = 0.2; p^ = 0.7. This is the case so far 
considered in literature [3T1I28]. For this case, by taking d 1 — 02 V "^1 
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This way, a density ratio of 1 : 1000 can be obtained, that the asymptotic behaviour is controUcd just by the 

as shown in figure [2l where a droplet of fluid of density ratio /ii//i2 and is independent oi pi, p2- Thus, keeping a 

P2 = 1 is in equihbrium with a surrounding second com- density ratio of 3, we can describe different binary fluids 

ponent of density pi — 1/1000. Notably, this conflguration changing appropriately the value of kinematic viscosity in 

has been obtained with ai = 0, a2 = 0.998 with a strength order to get the correct ratio between the dynamical vis- 

of the interaction potential imposed to be 512 = 321 = 0.5. cosities. This is sufficient to obtain a correct macroscopic 

Unfortunately, this is valid in a static case, but it is not description of the flow. In particular, we have chosen: 
stable for dynamic ones, since the corresponding value of 



pi = 0.34; p2 = 1-; vi = 0.042 V2 = 0.667 (18) 

and, as described above. 



c^ is so low that this component is affected by high-Mach 
number effects and, thus, cannot be correctly described 
in terms of our LB model. Furthermore, we have carried ^^^ — 4/9- q,^ — 0.8' 012 — 021 = 0.065. (19) 



out numerical simulations varying the density ratio and 
we have found that the highest density ratio attainable in 
our framework without the arising of instabilities is 1 : 3. 
This result is in agreement with what was found in a re- 
cent article devoted to this issue j36], but it is the first 
time that a complete dynamical test-case is considered. 
More specifically, the values chosen for this configuration 
are: ai = 4/9, a2 = 0.8 with the strength of interaction 
potential imposed to be 512 = (721 = 0.065. Recently, it 
has been pointed out that higher density ratios may be 
obtained also in microflows considering self- and cross- 
collisions [37]. Nevertheless, this strategy is based upon 
complex modelling and interpolation techniques, which 
could result computationally too demanding for general 
dynamical cases. Therefore, in order to model the usual 
situation of two fluids with very different densities as air 
and water or water in equilibrium with its vapour (in these 
cases the density ratio is « 500), we propose to work on 
viscosity values. Looking at eq. P^ . it is readily seen 



This choice corresponds to a dynamical viscosity ratio 
which can approximately describe the case of the capillary 
filling of water in air. In this conflguration the surface ten- 
sion has been computed to be 7 w 0.02 via the Laplace's 
law which states that at equilibrium 7 — -# for a drop of 
radius R. 

We have seen that the desired wettability, that is a 
given contact angle, can be introduced in our model through 
a phenomenological force [3T1I38] , whose expression is given 
by (fTT|) . Thus, the interaction with the walls is charac- 
terised by a phenomenological force proportional to giwRw 
and g2wPnn respectively for component 1 and 2, where giw 
and g2w are independent. Therefore, two different scenar- 
ios are possible for a given contact angle. Referring to 
fig- n we can let interact the incoming fluid 1, by im- 
posing a non-zero value of giwRw and giving to the force 
pT|) the sign to make the force attractive. On the other 
hand, it is also possible to impose a non-zero value of 
g2wPw and therefore to let interact the wall with the com- 
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ponent 2, but with the opposite sign, such that in this urc[5l There is a good agreement between analytical and 
case the fluid 2 is repelled. In both cases, the fluid 1 numerical solutions. 

feels hydrophilic walls with respect to the other compo- 

Finally, we have investigated the existence and the be- 

nent. However, in one case the fluid I is directly pushed 

haviour of precursor films, whose presence is indicated by 

inside by the attraction of walls, in the other it enters, 

a thin layer of liquid near-to-the-wall which penetrates the 

since the walls are hydrophobic with the fluid 2. In fig. 

capillary ahead the bulk meniscus. The presence and the 

m we show that for the same geometric and physical con- 
importance of precursors film has been first recognised by 

figuration and with an appropriate value of gkw and pw 

De Gennes [I] and much attention has been devoted to it 
the same behaviour is found and therefore from a macro- 

[40] . especially for the case of a spreading drop. Some ex- 
scopic point of view these two methods can be considered 

periments have been attempted even at nanometric scale 
equivalent. In particular, for this experiment we have used: 

[41j . In our numerical simulations, we have observed that 

giw = —0.065; g2w — 0; and p^ — 0.5; for the first case 

a thin layer (1 — 2 grid points) of liquid attached to the 
and giw — 0; g2w — 0.065; p^ — 0.7; in the second. This 

walls tends to penetrate faster, indicating the formation of 

choice corresponds to a contact angle 9 ^ 35. Data show 

a precursor film. The density of the film tends to decrease 

a good agreement with the asymptotic law given by the 

with the distance from the bulk of the liquid, since only 

Washburn's law, eq. ([T]), computed with our physical pa- 
some "molecules" escape from the bulk and they can not 

rameters and this contact angle. 

create a film as dense as the bulk liquid far away form the 

Then, we have concentrated our attention on the case bulk interface. We define the end of this film at the point 
of complete wetting. In order to obtain this configuration where the density of the liquid becomes 1/3 the bulk den- 
it is necessary to strengthen the interaction with the wall. sity. In figure [6l a snapshot of the density is shown which 
In figureini we show results for the different configurations supports visually the presence of this precursor film. It 
treated. For the values of physical parameters described is interesting to note that the layer is almost limited to 
in eq. ([T8|) . we found that the value of the contact angle the first point near to the wall, as observed in molecu- 
increases with increasing p^ until the value of pw — 1.0. lar dynamics simulations [29j . In figure [71 we analyse the 
After this threshold is reached, all profiles collapse. In fig- dynamical behaviour of precursors. In figure [71 the front 
ure we show the curves corresponding to p^j = 0.7, 1.0 displacement in time of both the front at the centre of the 
and pu, ~ 1.2. This behaviour indicates that the asymp- interface between the two components and the precursor 
totic value 6* = is reached. Then, we have computed the (as defined above) are shown. It is clear that the precursor 
solution of the equation (|15p for the parameters chosen film advances ahead the interface with a different law. Af- 
and = 0. The resulting solution is also shown in fig- ter 30000 time-steps the precursor has attained the length 
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of about 20 Ax. This behaviour has been further analysed nents and different kinematic viscosities have been used, 

and the difference between the precursor position and the such that the ratio between dynamical viscosities is com- 

interface one is displayed in figure [51 Both fronts advance parable with that for a water-air mixture. Furthermore, 

in time according to a square-root law, but with different it has been shown theoreticaUy that this ratio is rcspon- 

pre-factors. The displacement of the difference has been sible for the dynamical behaviour of the filling, fn this 

fitted with a square- root with a pre- factor of a = O.f . In sense, we have studied two important asymptotic config- 

the very early stage (t < 5000), the numerical curve of the urations: one one hand, when the two components have 

difference shows a kink indicating that it advances faster equal density and viscosity, the front advances linearly in 

than according to a square-root law, see fig. [H This can time. On the other hand, for a vanishing ratio fJ.2/ fJ'i, the 

be probably explained by the fact that, during this initial washburn's law is retrieved, which states that the front 

time, the bulk front is affected by a "vena contracta" ef- advances with a square-root law in time. Both asymptotic 

feet, which reflects the non-trivial matching between the behaviours can be found out through the analytical solu- 

reservoir and the capillary dynamics at the inlet [42]. tion of the equation of motion eq. p^ . 

These results appear to be in line with recent molec- j^ a configuration of complete wetting {9 = 0) the 

ular dynamics findings [29] and also with experiments at displacement of the interface with time has been shown to 

macroscopic scales [43] . agree with theoretical expectations. 

Finally, the presence of a precursor film ahead the cen- 
tral front has been detected. This film occupies only a thin 
In this work, we have studied the penetration of one fluid layer of one-two grid points and it moves with a square- 
into a capillary initially filled by another immiscible fluid ""^^^ t™e law (like the interface at the centre of the chan- 
by the means of a Lattice-Boltzmann model for binary ^el), but with a different pre-factor, faster than the inter- 
fluids based upon a pseudo-potential interaction. face at the centre. These resuhs appear to be in fine with 

A phenomenological force has been introduced in order 'recent molecular dynamics findings, 

to let the interface to develop a contact angle with the In this work, we present results for a chosen 2-D chan- 

walls. Two equivalent ways of imposing a given contact nel configuration. Even though recent 3-D molecular dy- 

angle have been analysed. namics results seem to suggest that the landscape is not 

In order to reproduce the classical experimental con- much affected by the change of geometry, it will be im- 

figuration, where a much denser fluid penetrates into a portant to simulate the imbibition in different and more 

much lighter fluid (vacuum or air), the model has been complex geometries in order to study the dependence of 

modified to let a density ratio between the two compo- the phenomenon on the simulation parameters. In partic- 
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Fig. 1. Geometrical set-up of the numerical LBE. The two di- 
mensional geometry, with length 2L and width H , is divided in 
two parts. The left part has top and bottom periodic bound- 
ary conditions such as to support a perfectly flat gas-liquid 
interface, mimicking a "infinite reservoir" . In the right half, of 
length L, there is the true capillary: the top and bottom bound- 
ary conditions are those of a solid wall, with a given contact 
angle Q. Periodic boundary conditions are also imposed at the 
west and east sides. The figure is based on a figure in [27] . 

ular, the behaviour of the precursor films could be more 
complex. 
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Fig. 2. Single-droplet total density profile. The density ratio 
between the droplet and the surrounding fluid is 1000, while 
the interface remains limited to a few grid points. 
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Fig. 8. The displacement of the precursor with respect to the 
central interface is plotted versus time. The curve is fitted with 
a square-root law at^''^ with the coefficient that turns out to 
be a = 0.1. This means that both fronts, the central and the 
precursor, advance in time with a square-root law but with 
a difi'erent pre-factor. The precursor advances in time more 
rapidly. 



